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Abstract. The 3-dimensional coherence matrix is interpreted by emphasising its 
• invariance with respect to spatial rotations. Under these transformations, it naturally 

^ , decomposes into a real symmetric positive definite matrix, interpreted as the moment 

of inertia of the ensemble (and the corresponding ellipsoid), and a real axial vector, 
'■ corresponding to the mean angular momentum of the ensemble. This vector and tensor 

^ I are related by several inequalities, and the interpretation is compared to those in which 

O ' unitary invariants of the coherence matrix are studied. 

Q^' PACS numbers: 42.25Ja, 42.25Kb, 02.50-r 



1. Introduction 



In the standard theory of partial polarization in paraxial light ^1121, the 2x2 hermitian 
0\ ■ coherence matrix (with unit trace) is decomposed into components with respect to the 
^ ■ Pauli matrices. These components, the Stokes parameters, summarise the second order 
O . statistical information about the ensemble; in particular, the sum of their squares is 1 
^ ■ for a pure polarization state, and for a completely unpolarized ensemble. 
' ^ . There has recently been a revival of interest in the corresponding coherence matrix 

in nonparaxial light, where in general there is no well-defined propagation direction, and 
the hermitian coherence matrix is3x3|HlllllSlinilZlElin]- In these treatments, by analogy 
with the 2-dimensional case, generalised Stokes parameters are defined by decomposing 
^ ■ the coherence matrix with respect to the Gell-Mann matrices; a generalized degree of 
. polarization jHl IH El El EI may be defined using the sum of squares of these components. 

Here, I propose a complementary interpretation of the 3x3 coherence matrix, 
motivated by geometric reasoning. Rotational, rather than unitary, invariants of the 
coherence are emphasised, and the matrix is found to decompose into its real part, 
which is symmetric and interpreted geometrically as an ellipsoid, and its imaginary 
part, which is antisymmetric and equivalent to an axial vector. The ellipsoid and vector 
have natural interpretation in terms of the ensemble of polarization states, and are 
related by certain inequalities to be described. 

Pure states are represented by a complex vector E, representing the electric field, 
in either two or three dimensions. This is represented geometrically by an ellipse by 
taking Re{-Eexp(— 1%)} and varying x (this may represent time evolution) [6| ITn| ITT|: 
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the ellipse therefore has a sense of rotation. In two dimensions, this is taken in the 
natural sense with respect to the plane, and polarization is either right- or left-handed. 
In three dimensions, the plane of the ellipse may vary, and the sense of rotation is 
a direction normal to the ellipse, defined in a right-handed sense with respect to the 
ellipse rotation jTHl HHj- The eccentricity of the ellipse can be 1 (corresponding to 
linear polarization), (corresponding to circular polarization), or any value in between. 
The ellipses are normalised in units of intensity \E\^. Polarization ensembles may be 
visualised geometrically as the set of polarization elhpses in the ensemble, adding 
incoherently. 

The paper proceeds as follows: the following section is a review of conventional 
2-dimensional coherence matrix theory; in section El the geometric decomposition of 
the 3x3 coherence matrix is described; section HI is devoted to the properties of the 
coherence matrix, and sectional to examples for certain ensembles. The paper concludes 
with a discussion in section IHl 

Polarization coherence matrices are special (classical) occurrences of density 
matrices, perhaps more familiar in quantum mechanics ^21 (pure polarization states 
corresponding to pure states, etc). Standard properties of density matrices (i.e. positive 
definite matrices with unit trace) will be employed without proof. 



2. The two-dimensional coherence matrix 



This section is included as a comparison for the 3x3 case, and reviews standard material 
discussed, for example, in [HEHEl- 

The 2-dimensional coherence matrix p2, assumed normalised (i.e. trp2 = 1), is 
defined 

[ {EyED {EyED ; ' 

where (•) denotes ensemble averaging over the ensemble of 2-dimensional complex 
vectors E = {Ex,Ey). p2 is normally expressed in terms of the Stokes parameters 
5*1, 5*2, 5*3, which are the components of p2 with respect to the Pauli matrices: 

The three Stokes parameters may be written as a 3-vector, the Stokes vector 

P= {81,82,83) (2.3) 

whose length |-P| is written P. 

Being a density matrix, p2 is positive definite (its eigenvalues are nonnegative), so 

det P2 = (1 - 8f - ^2' - 5|)/4 > 0, (2.4) 

that is, 

P < 1, (2.5) 



Geometric interpretation of the 3 -dimensional coherence matrix 



3 



which geometrically restricts P to lie within a sphere of radius 1, the Poincare 
sphere. This fundamental inequality is more commonly derived using the equivalent 
fact trp2 < (trp2)^- 

If the ensemble represents a single state of polarization (i.e. a 'pure state'), the 
coherence matrix is idempotent, 



Taking the trace implies that Ppurc = 1- On the other hand, if the ensemble is completely 
unpolarized, so p2,un is 1/2 times the identity matrix, then Pun = 0. This leads to the 
important decomposition of p2 into pure and unpolarized parts. 



The state of polarization corresponding to p2,pure here is the eigenvector corresponding 
to the larger eigenvalue of p2, and 1 — P is twice the smaller eigenvalue. 

The previous statements justify P as the degree of polarization. It is invariant with 
respect to any unitary transformation up2u''', by ()2.4|) (here, u represents an arbitrary 
2x2 unitary matrix). By the well-known relation between 2x2 unitary and 3x3 
orthogonal matrices, such unitary transformations correspond to rotations of the Stokes 
vector P. The operation of a unitary transformation on polarization states (or their 
ensemble average) is physically interpreted as the operation of a phase retarder [Hj, and 
the degree of polarization is unchanged when the ensemble is passed through a retarder, 
or series of them. 

The Stokes vector ()2.3j] resides in an abstract, 3-dimensional (Stokes) space, and the 
representation of phase retarders by 3-dimensional rotations is correspondingly abstract. 
If p2 is transformed by 2-dimensional rotations, corresponding to a real rotation of the 
transverse plane (i.e. op20^, with o 2 x 2 orthogonal), 5*1 and 5*2 may change keeping 
5*^ + 5*1 constant; S3 remains unchanged. An example case is the rotation in which Rep2 
is diagonalised: 



For pure states, for which the Poincare sphere representation is useful, the Stokes 
parameters provide geometric information about the polarization ellipse [HIE]- Si and 
5*2 inform about the alignment of the ellipse axes, the major axis making an angle 
arg(5'i + 182) /2 with the x-axis. S3 gives the ellipse area 1183, signed according to 
polarization handedness, so 5*3 is zero for linear, and ±1 for circular polarization. 
Obviously, 2-dimensional rotations only affect Si and 5*2; the rotation giving ()2.8j] 
represents aligning the major ellipse axis along x, the minor along y. 



P2,pure ~ P2,pure- 





(2.7) 
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3. Geometry of the 3-dimensional coherence matrix 



The 3-dimensional coherence matrix p = p3 is analogous to ^2.1^ . but with E = 

/ {e.,e:) {e,e;) {e^e:) \ 

p= {EyED {EyE;) {EyED . (3.1) 

\^ {E^ED {e^e;) {e^ed ) 

As before, it is assumed that trp = 1. 

The Gell-Mann matrices ^3] are the generators of 3-dimensional unitary matrices, 
just as the Pauli matrices generate 2-dimensional unitary matrices. Therefore, the 
generalised Stokes parameters Aj, i = l,...,8 [1111113111111, may be defined 

1 + A3 + Ag/VS Ai - iAs A4 - iAs \ 

Ai + iAs l-Aa + Ag/Vs Ag - iAy . (3.2) 
A4 + iAs Ae + iA7 1 - 2A8/V3 j 

(Other accounts, such as jH], use a different set of generators.) The analogies between 
fl2.2|l and ()3.2|1 are obvious: A3 and Ag, only appearing on the diagonal, generalise 5*1; 
the terms in the symmetric, off-diagonal part, Ai,A4,A7, generalise 5*2; and A2,A5, A7, 
appearing in the antisymmetric, imaginary part, 6*3. In particular, if A4,...,A7 = 
and Ag = then the remaining parameters are proportional to the usual Stokes 

parameters. This motivates the definition of the generalised degree of polarization P3 
HIHlin] as 



^=3 



(3.3) 



(A slightly different form was defined by |3](||5j.) This is the length of the 8-dimensional 
generalised Stokes vector, and, just as in the 2-dimensional case, it is invariant with 
respect to 3 x 3 unitary transformations. Since (trp)^ — trp^ > 0, it is readily shown 
that < P3 < 1. 

Although this approach is mathematically correct, it is not clear physically what 
P3 represents. Unlike the 2x2 case, in which the Stokes vector represents the complete 
state of polarization using three dimensions (which is easily visualised), the generalised 
Stokes vector requires eight dimensions, which is not so intuitive. 

There is a more serious problem with treating the 3x3 coherence matrix completely 
in analogy with the 2x2 case - there is no obvious physical interpretation via optical 
elements of 3 x 3 unitary transformations (nor any corresponding nonparaxial Jones 
or Mueller calculus). In two dimensions, as an ensemble of plane waves with the 
same direction but different polarizations propagates through an optical element, the 
corresponding coherence matrix is transformed by the appropriate Jones matrix, which 
is unitary for a retarder. In three dimensions, the ensemble of plane waves averaging to 
the 3x3 coherence matrix do not share a common propagation direction in general; any 
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physical device, represented by a 3 x 3 unitary transformation, should be insensitive to 
the propagation directions of the separate members of the ensemble. Mathematically, 
it is possible to find a unitary transformation which takes any 3-dimensional state of 
polarization E = {E^, Ey, Ez) to any other (leaving constant); there is no obvious 
physical situation in which different states of polarization in three dimensions undergo 
the same unitary transformation. 

It is physically and geometrically natural, however, to consider under orthogonal 
transformations rather than unitary ones; if viewed as passive rotations, this is simply 
equivalent to redefining cartesian axes in 3-dimensional space, and no physical operation 
at all. Clearly, under rotation, where p becomes opor^ (o 3 x 3 orthogonal), the real 
and imaginary parts of p transform independently of each other. The real part is a 
positive definite symmetric matrix with five parameters Ai, A3, A4, Ag, Ag. Since the 
(unit) trace is also unaffected by rotation, it may be considered as distinct from the 
rest of the real part. The imaginary part is a real antisymmetric matrix with three 
parameters A2,A5,A7, and in fact the triple (A7, — A5,A2) transforms under rotation 
like an axial vector (noted in [7j). p therefore decomposes into three parts: a real scalar 
(the trace), a real axial vector, and a real traceless symmetric matrix. These different 
parts (scalar, vector, tensor) are called irreducible tensor operators in group theory; 
the same decomposition occurs for density matrices of atoms with quantum spin 1, for 
which the vector part is called the orientation, the tensor part the alignment jTHj. 

From an analytical viewpoint, it is convenient to represent p using cartesian axes 
xi,X2,X3 with respect to which the tensor part is diagonal, giving 

/ Ml -iNs m2 \ 



P 



(3.4) 



\ -iN2 iNi M3 I 

The diagonal elements of ()3.4|) are restricted: 

Ml + M2 + M3 = 1, 1 > Ml > M2 > M3 > 0, (3.5) 

which follows from the fact that the tensor M = Rep is positive definite. It is 
geometrically convenient not to separate the scalar and (traceless) tensor parts of p, 
and this is not done in ()3.4|1 . ()3.4j) is analogous to ()2.8|) : the real part M has been 
(passively) diagonalised, leaving an off-diagonal imaginary part, which transforms as an 
axial vector 

N={Ni,N2,N^) (3.6) 

{\N\ is invariant under rotations). M and N have a simple geometrical interpretation, 
as follows. 

The real symmetric matrix M may be interpreted as the moment of inertia tensor of 
the ensemble. Geometrically, it is the moment of inertia of the set of polarization ellipses 
in the ensemble (taking each as an elliptical ring with uniform mass per unit length. 
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insensitive to the ellipse handedness) . As with moment of inertia tensors in mechanics, 
it may be represented in terms of its inertia ellipsoid, whose points {xi,X2-,x^) satisfy 

,y>2 /-v»2 /v»2 

The ellipsoid axes are aligned in the 1,2,3 directions, with lengths -^/Mi, ^/M^, VM^. If 
M3 = 0, the ellipsoid is flat (a;3 = 0). In general, the inertia ellipsoid is specified by 
6 parameters (the trace and Ai, A3, A4, Ag, Ag); the diagonal form in ()3.4|) . with three 
parameters, reflects that three Euler angles have been used implicitly in the choice of 
axes 1,2,3. The traceless part, dependent on the A parameters only, gives a measure of 
departure of this inertia tensor from isotropy. 

The vector N also has a simple interpretation as half the expectation value for 
(spin) angular momentum in the ensemble, 

5 = tr(Sp) = 2Ar, (3.8) 

where the spin matrices Sj for spin 1 in a cartesian basis are given componentwise by 
Si,jk = —i^ijk [m Hn], with Eijk the antisymmetric symbol. The axial vector N is 
therefore an average of the angular momentum, that is, the average sense of rotation of 
the ellipses, in the ensemble. Its direction, in general, has no relation to the principal 
axes of M (although its maximum length is limited by them, as described in the next 
section) . 

The inertia tensor M and orientation vector N therefore provide information about 
the real, 3-dimensional geometry of the polarization ensemble, and they rotate rigidly. 
Under more general unitary transformations (which have no physical interpretation), 
the eigenvalues of M and components of N may change arbitrarily (although keeping 
the unitary invariants tr p, tr and det p fixed). 

As an example of the geometric interpretation, figure Q is a representation of the 
inertia ellipsoid, orientation vector and dual ellipsoid (defined in the next section) for 
the matrix 

/ 14 -2i 2i \ 

(3.9) 



1 



Pex 20 



2i 5 -i 
y -2i i 1 y 

p may also be represented by its eigenvectors; if pa, pb, Pc represent the pure, 
idempotent coherence matrices corresponding to the eigenvectors of p with eigenvalues 
Aa, Ab, Ac (i.e. the principal idempotents 0^), then 

P = Kpa + \pb + AcPc- (3.10) 

The 2-dimensional analogue to (j3.1(jp immediately gives rise to the decomposition (j2.7p . 
Since the decomposition in ()3.10|) is in terms of three density matrices, a decomposition 
in terms of a single purely polarized part and unpolarized part is, in general, impossible 
(as previously noted in lUE])- The eigenvalues of p, being unitary invariants, do not have 
a simple geometric interpretation in terms of the inertia ellipsoid or orientation vector; 
however, the set of three eigenvectors rotates rigidly. This eigenvector representation 
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Figure 1. The inertia ellipsoid (black mesh), dual ellipsoid (grey mesh) and orientation 
vector corresponding to Pcx; in the xi^X2,x^ frame. Here, the orientation vector lies 
inside the dual ellipsoid, and not on its surface. 

of p provides a different geometric representation to that given by the inertia elhpsoid 
and orientation vector. However, it not not geometrically obvious when a given triple 
of polarization ellipses in three dimensions represent orthogonal polarization states, 
and solution of cubic equations is required to find the eigenvectors; moreover, the 
eigenvectors are not uniquely defined at a degeneracy. By comparison, the inertia tensor 
and orientation vector may be extracted directly from p and are always unambiguously 
defined. 

4. Inequalities satisfied by p 

In this section, various inequalities for M and N shall be found, using the fact that p 
is a statistical density matrix. 

Firstly, the Cauchy-Schwartz inequality may be applied to the off-diagonal elements 
of p in giving expressions of the form 

\{E^E;)\'<{\E^\^){\Ey\^). (4.1) 

Using the representation ()H.4j) . these imply 

Nl < M2M3, Nl < M1M3, Nl < M1M2. (4.2) 

Geometrically, this implies that the orientation vector N is confined to a cuboid with 
vertices {±^/M^,±^/M^,±^/M^). Since p is a density matrix, (trp)2-trp2 > 0, 
that is 

+ Nl + Nl < M2M3 + M1M3 + M1M2, (4.3) 

which is the sum of the inequalities ()4.2|1 . and therefore is less strong, geometrically 
restricting N to lie within the sphere circumscribing the cuboid defined above, 
(trp)^ — trp^ > 0, which is the distance by which N fails to touch the surface of this 
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sphere, is a unitary invariant. The traces of higher powers of p satisfy other inequahties, 
such as tr < tr p^ tr p, but such inequahties can be shown to be consequences of ()4.2|) . 
Nonnegativity of det p imphes that 

MiNl + M2NI + M^Nl < M1M2M3. (4.4) 

If M3 ^ 0, then 

Nl Nl A^l , 

\ \ — < 1 (4 5) 

M2M3 M1M3 M1M2 ~ ' ^ ' ' 

which geometrically means that N lies within the ellipsoid with axes in the 

1,2,3 directions, and lengths VM2M3, VM1M3, VM1M2. This ellipsoid is therefore 

circumscribed by the cuboid ()4.2|1 . and ()4.5j) is a stronger inequality than ()4.2|1 . The 

relationship between this ellipsoid and the inertia ellipsoid (|3.7|) justifies calling this 

ellipsoid the dual ellipsoid. If M3 = 0, (jOjl implies that A^i = iVg = 0, and (ji?^ gives 

|A^3| < A/M1M2; if the inertia ellipsoid is fiat, the dual ellipsoid is a line normal to it. If 

Ml = 1, M2 = M3 = 0, then the inertia ellipsoid is a line and N = 0. 

As with (I2.4jl . the fundamental inequality for the 3x3 coherence matrix is 

nonnegativity of the determinant, which is stronger than inequalities constructed using 

the trace. The geometric interpretation of the unitary invariant det p is the product 

of the distance by which N fails to touch the dual ellipsoid with the dual ellipsoid 

volume. This quantity, the trace, and the invariant discussed above are the only unitary 

invariants of p. Unlike the 2-dimensional case, the properties of p are complicated by 

the fact that polarization information is contained within both the inertia ellipsoid M 

and the orientation vector N. 



5. Examples of 3 x 3 polarization ensembles 



Completely unpolarized waves in three dimensions are a common occurrence, for 
example black body radiation. In this situation, the 3x3 coherence matrix is the 
completely unpolarized matrix pun, equal to one third times the 3x3 identity matrix 
(and P3 = 0). 

Coherence matrices for pure states of polarization satisfy Pp^-e ~ Ppure- 

Using (Q 

and ()3.5j) . this imphes that 

/ Ml 

Ppure 



v 







-iiV3 
M2 





0/ 



(5.1) 



with |A^3| = \/MiM2, M1 + M2 = 1. This is equivalent to a pure state in two dimensions, 
and represents a polarization ellipse E = {^/M^,±i^/M^,0) in 1,2,3 coordinates. The 
ellipse major axis is in the 1-direction, the minor in the 2-direction, and N is normal 
to the plane of the ellipse (oriented in a right-handed sense of rotation around the 
ellipse). If Ml = M2 = 1/2 in (jS.lj) . the state is circularly polarized, and N3 = ±1/2. If 
Ml = 1, M2 = (implying A^3 = 0), it is linearly polarized, det Ppure is zero, but unlike 
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_ 1 

Pexl — 2 



the 2x2 case this is not a sufficient condition for a pure state in generah trp^ must 
also be 1. The inertia elhpsoid of (jS-ip is flat, and N lies on the 'surface' of the (linear) 
dual ellipsoid, with equality in ()4.2j) . 

If the state is not pure but M3 = 0, then p satisfies fjS-H) with Mi + M2 = 1, but 
lA^sl < \^MiM2. An example is the density matrix 

/ 2 \ 

10 . (5.2) 
^ y 

The inertia ellipsoid here is fiat, and A?" = 0. It cannot be a pure state since 
Pexl 7^ Pexl • This matrix provides an example of a 3 x 3 coherence matrix which cannot be 
decomposed into the sum of a pure polarization matrix and the completely unpolarized 
matrix, since there is a zero on the diagonal - pexi — apun, for any positive a, leaves a 
matrix which is not positive definite. 

It is easy to visualise ensembles which have N = 0: their average angular 
momentum is zero. This may be achieved, for instance, by requiring for every E in 
the ensemble, E* has the same statistical weight as E. pexi is therefore the coherence 
matrix for the ensemble consisting of the pair of states (with equal weight) 

^cxi = {(V2,i,0),(V2,-i,0)} (5.3) 

(of course, this ensemble is not unique in averaging to pexi)- The ellipses corresponding 
to the pair ()5.3p are identical apart from their senses of rotation, which are opposite. 

Pexl is an example of a coherence matrix with N = 0, although its M is not 
isotropic; that is, the shape of the inertia ellipsoid is not constrained by the direction 
of the orientation vector. More surprising, perhaps, is that the converse is true - the 
inertia ellipse may be isotropic yet N takes on the maximum value allowed by ()4.2p . for 
example 

1 -i \ 

(5.4) 

/ 

which is the sum of pun and a completely antisymmetric matrix (which is not a density 
matrix). An ensemble which corresponds to Pex2 is the pair of states with equal weight 

£:e.2 = {(l,i,-l),(l,i,l)}. (5.5) 

The ellipses represented here share their minor axis (in the ^/-direction) and have 
orthogonal major axes. They both have the same shape (eccentricity 1/^/2), which 
geometrically implies that their total moment of inertia is isotropic (higher averages 
than quadratic are not isotropic). This pair of ellipses, along with the spherical inertia 
ellipsoid and orientation vector, are shown in figure El 

Both Pexl and Pex2 have the same eigenvalues 2/3,1/3,0 (equivalently, the same 
unitary invariants tr p, tr p^, det p); however, the two ensembles ()5.3|) . ()5.5|) are clearly 
not the same: the ellipses in the two ensembles have the same shape (eccentricity 




Pex2 = o I i 1 
1 
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Figure 2. The pair of polarization ellipses corresponding to the ensemble (15 .511 (grey), 
with their spherical inertia ellipsoid (black mesh) , spherical dual ellipsoid (grey mesh) 
and orientation vector, which here is vertical and on the surface of the dual ellipsoid. 



l/-\/2), but the orientations in space are different, and there is no obvious physical 
transformation between the two sets of states. 

In general, the minimum number of states in an ensemble required to specify p is 
three, and in fact, the (complex) eigenvectors of p suffice, as in fl3.10|l . In this case, 
the eigenvectors make up the ensemble, the probability weighting for each being the 
corresponding eigenvalue. Since Pexi,Pex2 each have one zero eigenvalue, an ensemble 
consisting only of two states is sufficient for these examples (the states in £^exi,i^ex2 are 
linear combinations of the eigenvectors, and are not orthogonal). 

6. Discussion 

Interfering nonparaxial polarization fields in three dimensions are more complicated 
than their paraxial counterparts, and their analysis involves subtle geometric reasoning 
[TUl ITHl IT^ ITT] . Most importantly, the Poincare sphere description breaks down for 
polarization states in three dimensions, because it cannot account for the direction of 
the ellipse normal iV; the appropriate nonparaxial analogue of the Poincare sphere is 
the Majorana sphere, which involves the symmetric product of two unit vectors, which 
describe the geometry of the nonparaxial polarization ellipse EUl ^j. These two 
vectors have a complicated expression in terms of the pure field state E. 

It would be of interest to find the relationship between the 3x3 coherence matrix 
and ensembles defined in terms of the Majorana sphere; a natural physical case would 
be when the Ex,Ey,Ez field components are gaussian distributed (for example black 
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body radiation). In this case, for a given p the distribution on the Majorana sphere 
would be unique and related to other gaussian Majorana statistics [21] • The analogous 
2x2 distributions on the surface of the Poincare sphere have a rather simple form 
[HI 122 123 1211 • Given the analytical complications of the Majorana sphere, it is unlikely 
that the 3x3 calculations will be straightforward, and it is unclear whether the geometric 
interpretation presented here would be helpful in this problem. 

As described in sectional there is no unique direction, or set of directions, associated 
with propagation for a general 3 x 3 p, and therefore, unlike the 2x2 case, there is 
no physical interpretation of 3 x 3 unitary transformations using conventional optical 
elements. Despite this lack of propagation information, a natural application of the 3- 
dimensional coherence matrix is in scattering theory, since a scatterer, such as a Rayleigh 
particle, responds only to the statistical E field at its position, i.e. the coherence matrix. 
It is therefore possible that classic problems such as atmospheric radiative transfer [22] 
may be analysed using the 3x3 coherence matrix. 

A natural experimental situation in which the nonparaxial coherence matrix is 
relevant is in the optical near field, for which measurements of the 3-dimensional field 
are possible [21! (^^ course the theory is not restricted to optical frequencies). The 
geometric interpretation should provide insight into the ensemble of polarization ellipses 
which gives rise to a measured 3x3 coherence matrix. 
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